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Abstract. This paper investigates quantum diffusion of matter waves in two-dimensional random poten- 
tials, focussing on expanding Bose-Einstein condensates in spatially correlated optical speckle potentials. 
Special care is taken to describe the effect of dephasing, finite system size, and an initial momentum dis- 
tribution. We derive general expressions for the interference-renormahzed diffusion constant, the disorder- 
averaged probability density distribution, the variance of the expanding atomic cloud, and the localized 
fraction of atoms. These quantities are studied in detail for the special case of an inverted-parabola momen- 
tum distribution as obtained from an expanding condensate in the Thomas-Fermi regime. Lastly, we derive 
quantitative criteria for the unambiguous observation of localization effects in a possible 2D experiment. 

PACS. 03.75.Kk Dynamic properties of condensates; collective and hydrodynamic excitations, superfluid 
flow - 42.25.Dd Wave propagation in random media - 72.15.Rn Localization effects (Anderson or weak 
locahzation) 



[ 1 Introduction 

■ In recent years an increasing number of theoretical and 
, experimental studies discussed the transport of ultra-cold 

■ atoms and Bose-Einstein condensates in the presence of 
disorder 1 (and references therein). One emblematic phe- 

. nomenon in the field is the celebrated Anderson local- 
' ization phenomenon, a disorder-induced metal-insulator 
, transition observed in the absence of inter-particle inter- 

■ actions |2|3|4j . Recently two papers have reported exper- 
imental evidence for exponential spatial localization of 
matter waves using dilute Bose-Einstein condensates in 
a ID optical speckle potential [5] and in a quasi-periodic 
optical lattice [6j. These results call for investigations in 
higher dimensions as it is known that transport in disor- 
dered ID and 2D potentials always occurs in the localized 
regime whereas in 3D there is a transition [7] . 

In two preceding articles '8"9' we presented a theoret- 
ical description for the quantum diffusion of monochro- 
matic non-interacting ultra-cold atoms in 2D and 3D dis- 
ordered optical potentials. In real experiments, a crucial 
question in this context is the influence of a finite initial 
momentum distribution of the atoms, together with the 
limits imposed by dephasing processes, like spontaneous 
emission, and boundaries. These questions are addressed 
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in the present paper, focussing on the case of 2D optical 
speckle potentials. Coherent diffusive transport of non- 
interacting matter waves is studied for any initial Wigner 
phase-space distribution within the framework of the self- 
consistent theory of localization [9110111] . Within this for- 
malism, we then analyze the influence of dephasing pro- 
cesses and of boundaries. The expected diffusion constant 
is calculated, both for the momentum distribution of an 
interaction-driven expansion of a Bose-Einstein conden- 
sate [12; and a Gaussian momentum distribution [13j . Fi- 
nally, we derive criteria for the experimental observation 
of Anderson localization in 2D correlated speckle poten- 
tials. 

The paper is structured as follows: we first recall in 
section[5]some basics about propagation in disordered sys- 
tems and calculate the expected probability density distri- 
bution and the variance of the expanding cloud. As a first 
application, we discuss the case of phase-incoherent, clas- 
sical diffusion in section [3l Section |4] considers coherent 
corrections to transport, viz., weak and strong localiza- 
tion. We calculate the diffusion constant self-consistently 
at finite frequency. In the stationary limit, this leads to 
the shape of the localized density distribution and the rel- 
evant localization length. In section [5l we incorporate the 
effect of dephasing by spontaneous emission into the for- 
malism, and derive criteria for the effective threshold sep- 
arating the diffusive from the (almost) localized regime. 
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where only very slow residual diffusion occurs. Sections [6] 
and [7] are devoted to the implications of finite system size 
and experimental realizability. 



2 Transport in disordered systems 

In this paper we describe the dynamics of non-interacting 
cold atomic gases evolving in static two-dimensional ran- 
dom potentials V{r). If interaction plays no role, a one- 
body description of the atomic dynamics is justified. The 
quantum dynamics is then described by the Hamiltonian 



"-t + ''''>- 



(1) 



where V{r) describes the potential fluctuations around its 
average value, which we choose to be the origin of ener- 
gies. In a single realization of the disorder, an initial plane 
wave k with kinetic energy Ek = h^k'^/{2m) will be ran- 
domly scattered by the potential fluctuations to all other 
accessible states. As a consequence, only statistical quan- 
tities obtained by configuration average over the disorder 
realizations can reveal generic transport properties. In the 
following, configuration averages will be denoted by (. . .), 
and as a first property of the disorder potential, V{r) = 0. 



The slow algebraic decay of the real-space correlations (as 
p~^) implies that their Fourier transform is non-analytic 
at K = 0: Popt('«) = 47r - 4|/«| -I- 0{k'^). Also the fact that 
the power spectrum has finite support in fc-space has im- 
portant consequences for the localization of matter waves 
in a ID geometry, allowing for the crossover from exponen- 
tial to effectively algebraic localization in position space 

EEH- 

Potential correlations are also often described (or ap- 
proximated) by a Gaussian correlator, such as in [TO], be- 
cause it is easily implemented in the numerics and leads 
to somewhat simpler analytical calculations: 

Cgauss(p) = exp (-/3^/2) , (5a) 

Pgauss('«) = (2^)^^/^ gxp (_k2/2) . (5b) 

If one wishes to match the small-r expansion of the 2D 
optical speckle (I3|) one must choose Cgauss = \/2Copt- The 
power spectrum (|5b|) decays as a Gaussian, which is often 
a qualitatively good approximation to the finite support 
of the true optical speckle spectrum (j4|). 

If only the behavior at very low momenta k = k(^ <ti 1 
is of interest, the detailed structure of the correlation func- 
tion is unimportant, and a good and simple approxima- 
tion is provided by P{k) w -P(O) where Popt(O) = 47r and 
Pgauss(O) — (27r)''/^. This in turn corresponds to an effec- 
tively 5-correlated position-space correlator 



2.1 Spatial potential correlations 



V{r + r')T/(r') = y^P(0)C^(r). 



(6) 



Experimentally, random potentials V(r) can be realized 
as optical speckle potentials generated by monochromatic 
illumination and imaging of an appropriate diffusing plate 
[14]. In the regime of parameters that we will explore, the 
pair correlation V{r + r')y(r') is the only ingredient one 
needs to calculate relevant microscopic quantities charac- 
terizing the atomic dynamics [2]. The potential fluctua- 
tions are then characterized by the pair correlation 



V{r + r')V{r'):^V'C{r/C) 



(2) 



that depends only on the modulus r = |r| when translation 
and rotation invariance is restored on average. The spa- 
tial correlation function C{p = r/Q) decays from C(0) — 1 
to zero over a characteristic spatial scale C,. In the case 
of an optical speckle potential this correlation length 
Copt = (afci)"^ is set by the wavelength Al = 27r/fcL of the 
monochromatic laser field and by the numerical aperture 
a of the imaging system (in current experiments a ~ 0.3). 
A 2D optical speckle potential created from a circular 
diffusive plate is described by [9115] 



Copt(p) - {'2Ji{p)/pf 



(3) 



where Ji {p) is the Bessel function of order 1 . Its 2D Fourier 
transform P(/«) = / dpe~^'^'''C{p), also known as the fluc- 
tuation power spectrum, is given by 



Popt(A^)=8 arccosf-f Jl-(f)2 0(2-k). (4) 



2.2 Weak disorder regime 

The correlation length Q defines a correlation energy 



mC^-^ 



(7) 



and its associated time scale r<^ = h/Ec^ = mC^'^ /h. The 
matter wave dynamics is then driven by three different 
energy scales: the kinetic energy Ek, the strength V of 
potential fluctuations, and the correlation energy Eq. 

The impact of potential fluctuations on the atomic dy- 
namics is simply estimated by the magnitude Vt /h of the 
random phase kick experienced by an atom travelling over 
the distance Q = vt &i velocity v = hk/m through po- 
tential fluctuations of strength V. As shown in [9 , the 
weak-scattering regime is realized if 






2Ek kQ 



Throughout the paper, we will use 



V 



V_ 
El 



(8) 



(9) 



as a measure for the potential fluctuation strength. Con- 
dition ([5]) shows that the weak-disorder regime is realized 
if the kinetic energy Ek is sufficiently above the important 
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energy scale Ea — V^/Eq. In 3D, this is essentially the 
mobility edge separating extended states with Ek > Ea 
from localized ones with Ek < Ea [9]- 

Even within the weak scattering regime E^ ^ Ea one 
may still be able to probe the low-energy or J-correlated 
potential regime Ek <^ E(^ provided the following energy 
hierarchy holds: Ea "^ Ek -^ E(^. This requires very weak 
potential fluctuations V <^ Ec^ ot equivalently ry <C 1. 

In the weak-scattering regime the atoms experience 
thus many random phase kicks of small amplitude in the 
course of time, which repeatedly alter their wave function. 
These phase kicks are the seeds of a diffusive behavior for 
the average probability density that will be calculated in 
the following. 



2.3 Phase-space dynamics 

The ensemble-averaged atomic dynamics in a d- 
dimensional potential is fully described by the average 
Wigner distribution [17 : 

l^t(r,k)= /dr'e-"^-"-' (r+^|^|r-^\ (10) 



where Q{t) = U{t) qq U'^{t) is the one-particle average den- 
sity operator. Here go denotes the initial atomic density 
operator and U{t) = 0{t) exp{—iHt/h} is the forward- 
time evolution operator for the Hamiltonian ([T]), being 
the Heaviside step function. From this Wigner function 
one can extract the marginals p(r, t) (spatial distribution 
at time t) and 7r(k, t) (momentum distribution at time t) 
according to: 

p{r,t) = J-^Wt{r,k), TT{k,t)= JdrWtir^k). 

(11) 
The normalization of these marginals is 

y^7r(k,t)=y'drp(r,t) = l, (12) 

which means that the Wigner distribution Wt is normal- 
ized as 

/|^H.-.(r,k).L (13, 

Rewriting the integrand of (jlOp by separating the ini- 
tial Wigner distribution Wq gives Wt in the form 

/dr'dk' 
^^-^Gt{r,k;r',k')Wo{r',k') (14) 

with a rather transparent physical meaning: starting from 
the point (r', k') at time t = 0, the initial quasi-probability 
Wq is propagated in phase space with Gt to give Wt- The 
phase-space propagation kernel Gt thus encapsulates the 
relevant quantum dynamics and reads: 

Gt^ J dridr2 e-'(''-i-k'--2) £^(r, r', n, rs) (15a) 

Ct = Xt(r+^,r' + f)if*(r-^,r'-f). (15b) 



Here the quantum propagator Kt{ri,r2) — {ri\U{t)\r2) is 
the probability amplitude to end up at point ri at time 
t when starting at point r2 at time i = 0. In Feynman's 
path description, Kt is expressed as a sum over all possible 
paths connecting r2 and ri. This assures that Gt indeed 
contains the non-local quantum interference that affects 
the density dynamics in phase space. It can be readily 
checked that Gt satisfies: 

lim Gt(r, k; r', k') = (271^ S(r - r') S{k - k') (16a) 



/i^^*('-'^'^''^') = ^' 



(16b) 



where the latter merely expresses, as it should, conserva- 
tion of the total density. 



2.4 Disorder-averaged probability density 

Using the previous results, it is easy to show that the 
average number density distribution is given by 

/dr'dk' 
.^^Ft{r-r'X}Wo{r',k'), (17) 

with a spatial propagation kernel given by 

/Air 
^^G,(r,k;r',k'). (18) 

Taking advantage of the translation invariance, one can 
develop the intensity kernel into its Fourier components 
at fixed momentum k, 

^cu, q, k) = f drdt e'('^*-'i'-) Ft(r, k). (19) 

Probability conservation and linear response theory [3] im- 
ply that for long times and distances, i.e., small lo and q, 
the kernel takes the diffusive form 



<?(w,q,k) 



1 



—iLU + D(llJ, k) q'^ 



(20) 



featuring the cj-dependent diffusion constant D{Lu,k) that 
can be calculated microscopically by quantum transport 
theory |9ll8j . In this paper, we only consider situations 
where the disorder average restores translational and ro- 
tational symmetry, as exemplified by the potential corre- 
lator 12]). This assumption implies in particular that the 
diffusion constant D(cj, k) must be a scalar and can only 
depend on |k| = k. Before turning to the microscopic the- 
ory, we first explore which arc the observable consequences 
of the diffusive form (|20p in transport experiments with 
cold atoms. 



2.5 Variance of the expanding cloud 

The position variance of the expanding cloud of cold atoms 
in the speckle potential gives direct access to the diffusion 
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constant. Indeed, the variance of the expanding cloud is 3 Classical transport: the Boltzmann picture 

defined by Ar^ = (r^) — (r) where (...) denotes the aver- 
age with respect to p(r,t). Using (fT7|) . (fTO)) . and (1^ one 
finds: 



Ar^{t) ^Arl^2d 



2^ 






(21) 



where Atq is the initial variance of the atomic cloud. 



The 



notation uj + iO indicates that the integration over fre- 
quency has to circumvent the pole at w = from above 
in order to describe forward-time propagation. Here, the 
effective diffusion constant 



V{u;)^{Diu;,k))^ 



dk 

(27rY 



7ro(k)L>(w,fc) (22) 



is the fc-dependent diffusion constant averaged over the 
initial momentum distribution 7ro(k). 



2.6 Interaction-driven momentum distribution 

We wish to calculate the average diffusion constant V^uj) 
using the momentum distribution of an expanding Bose- 
Einstein condensate (BEC) released from an isotropic har- 
monic trapping potential with frequency f2. Initially, the 
interacting condensate with chemical potential fi ^ hf2 is 
trapped in a Thomas-Fermi parabola with radius R de- 
termined by /i = ^mil^B?. When the harmonic trap is 
switched off, the expansion occurs in two steps. For small 
times I?i <C 1, the repulsive interactions drive the expan- 
sion and the condensate acquires a dynamical phase. For 
long times fit ':^ 1, the density drops so much that the 
expansion becomes free and the condensate wavefunction 
takes the form Hn] 



(r,t) ex vl — z^ 0(1 — z) exp(mz^). 



(23) 



with z — r / {\f2RQt) and u = 2fit/h. Since fit ^ 1 and 
IJ,/{hf2) » 1, one can use a stationary-phase approxima- 
tion with M ^ 1 to calculate the Fourier transform to the 
momentum domain, which is stationary: 



0(k)(xyr^(fclF6'(i-fcO- 



(24) 



Here, ^ = h/^Amfi is the healing length defined such that 
fi = h?/Am£^^. Hence, the normalized momentum distri- 
bution 7ro(k) = |0(k)p is given by 



7ro(k) 



(20'^ d{d + 2) 
~^d 2 



il~ikO^)0{l-kO (25) 



where Sd denotes the surface of the (d — l)-dimensional 
unit sphere (Si = 2,5*2 = 27r, ^a = An). This distribution 
has again the form of an inverted parabola with an upper 

cut-off A:,„ax = 1/^ M- ^ ^ 

When the asymptotic regime of free expansion is 
reached, the speckle potential is switched on. The distribu- 
tion ttq (k) then constitutes the initial momentum distribu- 
tion for the evolution of the condensate in the disordered 
potential. 



In order to put the above results to work, let us first 
calculate the dynamics of particles in a disordered po- 
tential if all interference effects can be neglected. This is 
the realm of classical transport theory, whose origins date 
back to the Drude-Boltzmann theory of metals, more than 
one century ago }18|20)21j . and that was later adapted 
to light propagation through interstellar atmospheres (ra- 
diative transfer theory) [22]. The basic physical ingredi- 
ents of this description are (i) first to assume that any 
possible interference effects are washed out under disor- 
der average (random-phase assumption) and (ii) second 
to devise a detailed-balance analysis of energy transfer 
in phase space (scattering, absorption, sources, etc). This 
powerful description leads to a physically very simple and 
appealing picture: sufficiently far from the boundaries and 
in the hydrodynamic regime of long times and large dis- 
tances, transport is described by diffusion. Note that the 
description of transport in terms of a Boltzmann equation 
as such — and by its diffusive limit at long time and large 
distances — is not restricted to the weak scattering limit 
established in section [^?^ above. However, it is in the weak- 
scattering regime that macroscopic quantities can be eas- 
ily computed from the microscopic ones as will be shown 
in the following. 



3.1 Boltzmann transport mean free path 

When all interference corrections to the intensity kernel 
can be neglected (see section [5] below for a more detailed 
derivation) one finds a finite diffusion constant in the sta- 
tionary limit of frequency w — > given by [S] 



Dn{k) 



md 



kiBik) 



viB 

d 



f2 
TB 



(26) 



Note that this expression formally remains valid beyond 
the weak scattering limit, provided the full dispersion re- 
lation (instead of the free-particle expression v — hk/m) 
is used for the average velocity. The Boltzmann transport 
mean-free path £b is the characteristic spatial scale beyond 
which memory of the initial direction is lost. It thus iden- 
tifies with the average step of the random walk induced 
by the scattering processes. This classical transport mean- 
free path can be calculated both from the microscopic the- 
ory [5] and from a detailed-balance Ansatz, by an angular 
integral over the differential cross section of a single scat- 
tering event, weighted by the scattering anisotropy factor 
(1-C0S6'): 



1 

ki^ 



rj- 



^4-d 



{2t:Y- 



y(1- COS 6l)P (2k sin(6'/2)). (27) 



Here k — k(^, df22 = dO in 2D (integration range from to 
2tt) and dfis — 2% sin 9 d9 in 3D (integration range from 
to tt). 

In the weak-scattering regime, the elastic scattering 
mean-free path £s is given by expression (|27|) with the 
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anisotropy factor replaced by 1. For the correlated poten- 
tials considered here (either the speckle or the Gaussian 
one), £b is larger than is, the two being asymptotically 
equal in the low-energy limit where scattering is isotropic. 
An analytic expression for £b is available in the small- 
and large-momentum limit [9] for the 2D optical speckle 
correlator ([3|) where C = Copt: 









457r(fcC)^ 
128772 



kC < 1, (28a) 
/cC>l- (28b) 



For the Gaussian correlation (j5bp with C — Cgauss, one 
finds in 2D an analytical result valid for all momenta fc, 



i:^-|fJe~'^^^^[^o(fcV)-/i(fcV)], (29) 



1 

Mb 



where I^ is the modified Bessel function. This result has 
also been derived in p3] in a different context. In the low- 
energy limit k( <^ 1 both optical speckle and Gaussian 
correlations result in the same expression 



1 



P(0)r/2 P{0)Ea 



k^C 



2Ek 



(30) 



just as substituting F (2Ksin(6'/2)) w P(0) in (gT]) would 
have given. Hence k£B S> 1 in the weak-disorder regime 
Ek > Ea- 

On the high-momentum side, it is important to note 
that contrary to the ID case [5|12j . in 2D the transport 
mean- free path £Bik) is finite for all fc- values, even though 
the potential correlator may have finite support, because 
even fast atoms can be deflected ever so slightly by a 
smooth random potential. 



3.2 The Diffuson kernel 

With the frequency-independent diffusion constant 
DB{k), the intensity propagation kernel is easily calcu- 
lated from (fTO|) and (|20p and reads 

FBir,k,t) = (47ri?B(fc)t)-'^/2 exp[-rV(4Z?B(fc)i)]. (31) 

This Boltzmann kernel, known as the diffuson [21j . obeys 
the diffusion equation: 



dtFB - DbV^Fb = S{r)S{t). 



(32) 



The diffuson obviously fulfills probability conservation. In- 
terestingly enough, it also boils down to 6{r) when i ^ 
even though its expression is in principle only valid at 
large enough times. The range of validity of this diffu- 
sive description is then expressed as r ^ ^b and t ^ tb 
in terms of the Boltzmann transport time tb = £b/Db- 
With the interaction-driven momentum distibution p4p . 
this will occur for the whole distribution as soon as 



Ant 
Hd 



>fcn 



>fcn 



(33) 



where £max = ^B(fcmax = I/O' 

Since diffusion in phase space occurs with a 
momentum-dependent diffusion constant, the average 
probability density p{r, t) itself does not obey a diffusion 
equation. But going back to (fT7)) . we note that for times 
large enough, the spatial width of the diffuson will always 
be much larger than the spatial width of Wo . This implies 
that, in the long-time limit, the probability density is well 
approximated by the momentum-distribution average 



p(r,i)«(PB(r,fc,i)), 



(34) 



Finally, with the frequency-independent diffusion con- 
stant DB{k), one finds that the variance pT|) of the ex- 
panding cloud increases as expected linearly in time. 



Ar^{t) ^ Arl + 2dVBt, 



(35) 



with the momentum-averaged diffusion constant analo- 
gous to (EH), 

Vb - {DB{k))^„ . (36) 

This average over the inverted-parabola distribution (f25|) . 
can be easily evaluated numerically for any of the mean- 
free paths (p7)l - p0)) above. An analytical result is available 
in the low-energy limit: 



Vr^ 



i?B(fcmax) _ 1 h 2^n 



3 2to PiO)EA 



(37) 



where i^max = 2/i in terms of the initial chemical poten- 
tial. This is a good approximation only if the momentum 
dispersion is so small that fcmaxC = C/C "^ 1j i-^-, in the 
limit where already the initial condensate is weakly inter- 
acting. 



4 Coherent transport and localization 

The ensemble average over static disorder alone cannot 
wipe out interference effects, and Boltzmann's description 
does not adequately describe phase-coherent propagation. 
Indeed, the constructive interference between amplitudes 
counter-propagating along loop- like paths survives the dis- 
order average and increases the particle's probability to 
return to its starting point. This increased tendency to 
stay behind translates into a reduced diffusion constant, 
an effect called weak localization (WL) |18I21I24] . Under 
suitable conditions, this can even completely inhibit diffu- 
sion, a phenomenon known as strong (or Anderson) local- 
ization (SL). It is rigorously proven that transport in ID 
bulk samples is always suppressed by strong localization, 
and on the basis of the single-parameter scaling theory 
this remains valid also in 2D, whereas in 3D localization 
occurs only when k£B ^ 1 (loffe-Regel criterion); further 
details can be found in the recommendable review [3]. 



4.1 "Loop-renormalized" difTusion constant 

Phase-coherent bulk transport in the weak localization 
regime fc^B S> 1 can be described by a diagrammatic 
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perturbation theory, developped by VoUhardt and Wolfle 
|10|18j . When weak- localization interference effects are 
self-consistently taken into account, the Boltzmann diffu- 
sion constant D^(k) is renormalized into an interference- 
reduced diffusion constant D{u)^ k) solving 



1 _ 1 Ad{k) 

D{u;,k) D^ik)^ Djiik) 



C{uj,k) 



(38) 



where Ad{k) = {27r)'^h/{mT:Sdk'^~^). This formula fea- 
tures the "return probabihty" C{uj,k), given by 



C{LU,k) 



dq 1 

(27r)'^ -iuj + D{LU,k)q'^ 



(39) 



in the diffusive regime. A mean-field hke, self-consistent 
description is obtained because the corrected diffusion 
constant D{uj,k) enters the diffusive propagator under the 
integral. The implicit equation for -D(w, k) can be conve- 
niently rewritten as 



Deik) ^ D{uj,k) 



ik<^- 



r.d-i 



dq 



q^ — iijj/D{ijj, k) 



(40) 



In our previous articles |8l9j . we have studied directly the 
stationary regime at w = 0. Presently, we will explore 
the consequences of the finite-frequency approach, which 
turns out to be quite powerful. 

The formal integral for the return probability (|40D may 
present divergencies, either infrared (for small q) or ultra- 
violet (for large q), that have to be regularized by suitable 
cutoffs. Working at finite frequency uj has the advantage 
that the integral stays finite as q ^ 0, so there is no need 
for introducing a special infrared cutoff. 

However, the return probability in 2D displays a log- 
arithmic ultraviolet divergence (g — > oo). We regularize 
it by introducing an appropriate upper bound (/max, de- 
fined as the minimal distance beyond which interference 
can play a role. Since the return probability is evaluated 
in the diffusive regime, we have chosen this minimal dis- 
tance to be the Boltzmann transport mean free path and 
accordingly imposed 1/^b as the upper bound. Please note 
that there is no need to dress also this ultraviolet cutoff 
by localization corrections: during the short time required 
to scatter through small loops, corresponding to oj k, t^ , 
the interference corrections will be shown to be negligible 
just below. 



4.2 Weak and strong localization in 2D 

Solving now (|40p for d = 2, we obtain the functional de- 
pendence oi X ~ D{co, k)/D'B,{k) on w at fixed k through 
the implicit equation |25j : 



X 



'lUJTB 



exp[7rHB(l-^)]-l' 



(41) 



In the weakly disordered regime fc^e ^ 1, it is advisable 
to rewrite this transcendental equation as 



X = 1- 



1 



irkiB 



log 



1 



X 



ILUTB 



(42) 




Fig. 1. Scattering mean free path £s, Boltzmann transport 
mean free path Eb [eq. (|27|l ] and localization length ^loc 
[eq. (144 |l ] as a function of the atomic wave vector k, all mea- 
sured in units of the correlation length (" = ^opt of a 2D speckle 
potential. The disorder strength is 77 = V/E( = 0.2. Note the 
extremely rapid increase of the localization length, a charac- 
teristic feature of 2D localization. 



which gives as perturbative solution the celebrated weak- 
localization correction |26l27j 



DB{k) 



= 1- 



1 



Trfc^T 



■log 



1- 



1 



IUJTb 



(43) 



In the frequency range lot-q ~ 1 of small diffusive loops, 
this correction is negligible. This means that the diffusion 
of a particle starts off with a classical Boltzmann random 
walk and is consistent with the choice of the unrenormal- 
ized ultraviolett cutoff gmax — 1/^b discussed in the pre- 
vious subsection. 

At intermediate times ujtb <^ 1 the weak localization 
corrections described by (jlSj) kick in. This expression is 
vahd as long as the correction is not too large, i.e., for 
ujtb > (-^b/Cioc)^ in terms of the new (and important) 
spatial scale 



60c (fc) = ^B(fc)Vexp[7rHB(fc)] - 1 « £Be""^/2 (44) 

whose significance will become clearer in a moment. For 
even smaller frequencies ojtb <C (fe/Cioc)^ the solution of 
(|4T|) crosses over to the asymptotic behaviour 



D{uj,k)^-tLo^l,ik). 



(45) 



This linear dependence on frequency with 
limtj_,o D{uj, k) = implies that the diffusion stops 
in the long-time limit. The intensity propagation kernel 
(fT9| calculated with the asymptotic solution (l45l) is 
stationary: 



Foo{r,k) 



iqr 



^q 

(2^)2 1 + q^^f 



Kn{r /£,\oc 



loc 



2<c 



(46) 



where Kq is a modified Bessel function 
totics for large distances Kq{x) 



The asymp- 
"^ shows an 



yjTT/2xe 

exponential decay, and we see that '^ioc(fc) is the localiza- 
tion length. 
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Fig. 2. Solid black line: stationary spatial probability density 
l|47p of a BEC wavepacket initially in the Thomas- Fermi regime 
at disorder strength 77 = 0.2 and fcmaxC = 0.65. Note the loga- 
rithmic scale and the fact that most of the wavepacket remains 
well localized near the initial position, an effect of the small-fe- 
components of the wavepacket. Only in the far tails where the 
density is extremely low one observes exponential localization 
with localization length ^ioc(fcmax) ~ 48 /im. The full shape 
of the wavepacket is well described by exp[— r/^ioc(fcmax)]/»'^ 
(red circles). 



As noted in section 13.11 above, in 2D the transport 
mean- free path £B{k) is finite for all fc-values. Conse- 
quently, the locahzation length ^\oc{k) as given by (|44p 
is also finite, and there is exponential localization for all 
fc-components in the initial momentum distribution ttq (k) . 
This implies that a wavepacket will follow a diffusive be- 
haviour at short time, but the dynamics will slow down 
and finally freeze. The resulting stationary probability dis- 
tribution depends on the initial state. Provided that the 
initial spatial size of the wavepacket is small compared to 
the effective width Ar{k) = ^ioc(fc) of the propagation ker- 
nel (for the k values populated in the initial wavepacket), 
the limit distribution reads 



lim p{r, t) 



Ko{r/^,Uk)) 



(47) 



This final density distribution is shown in fig.[2]for the typ- 
ical values 77 = 0.2 and fcmaxC — 0-65 used throughout this 
paper. Obviously, the shape near the center is not purely 
exponential. Only quite far in the tails is the wavepacket 
component at femax dominant, leading to exponential lo- 
calization with localization length ^ioc(fcmax) ~ 48 /im. 
Except very close to the origin, the density distribu- 
tion is quite well described by exp[— r/^ioc(A:niax)]/''^^^ 
as obtained by using the asymptotic behaviour of Kq in 
eq. (gZl). 

This behaviour in 2D can be compared to other di- 
mensions. In ID, and still within the Born approxi- 
mation, two regimes have to be distinguished. In case 
the highest momentum lies still inside the support of 
the speckle two-point correlator, all A:-components are 
exponentially localized, and an asymptotic behaviour 
exp(— |z|/^ioc(fcmax))/|-2r^^ is predicted [H]. In case the 
highest momenta are beyond the support of the two-point 
correlation, then a regime can be reached with an effective 
algebraic density decay as |z|~^; these predictions have 
been recently confirmed by experiments [5l. In 3D, aver- 
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Fig. 3. Solid curve: localization length ^loc [eq. (|44p ] as a func- 
tion of the atomic wave vector k for t] = 0.2, both in units of 
the 2D-speckle correlation length C,. Dashed line: stationary 
rms size Ar/2 of the wavepacket [eq. (|48|l ] resulting from the 
expansion of an initial Thomas-Fermi wavepacket with max- 
imum momentum fcmax = fc. As a result of the average, var- 
ious shorter localization lengths contribute and the resulting 
final size of the wavepacket is significantly smaller than in the 
monochromatic case with the corresponding k = fcmax. Note, 
however, that both curves display a similar rapid increase with 
fc. 



aging over the same momentum distribution leads to a 
|r|~^ algebraic decay, without any exponential envelope. 
This is due to atoms with momenta close to the mobility 
edge, at least for the limiting case of uncorrelated Gaus- 
sian disorder |29| . Thus, only in two dimensions do we 
have the favourable situation where the density profile al- 
ways shows an exponential localization. This is because 
the two dimensional localization length is always finite for 
all fc-components for correlated disorder even within the 
Born approximation. 

Solving next for the variance (I2ip we get: 



\iuiAv^{t) = Avl + 2d{^l,{k))^ 



(48) 



When the initial dispersion Ar^ can be neglected, the 
spatial variance will be given by the momentum-average 
of the squared localization length. This is the essence of 
strong localization: Transport starts in the Boltzmann 
regime and, as time increases, the diffusive dynamics gets 
slowed down by interference effects. Finally the dynam- 
ics freezes in the stationary limit and 2D transport al- 
ways reaches the localized regime within a bulk system in 
the absence of phase-breaking mechanisms. The time scale 
for the complete crossover from the diffusive to the sta- 
tionary localized distribution shown in Fig. [2] is given by 

the localization time Tioc(fcmax) = Cloc(^max)/-DB (fcmax) 

for the fastest atoms with momentum fcmax- Of course, 
some components of the atomic cloud (the ones with lower 
momenta) will localize earlier and on a shorter spatial 
scale, resulting in a strongly peaked probability density, 
see Fig. [2l 
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5 Quantum difFusion with limited phase 
coherence 

Since 2D localization relies on interfering amplitudes over 
loops of arbitrary size and since the bulk localization 
length ^loc increases exponentially fast with Mb-, one must 
imperatively discuss spatial or temporal factors limiting 
the interference effects. An important limitation is due to 
dephasing, and another obvious limitation is the finite size 
of the sample. Within the framework established so far, 
we are able to take into account both limitations. First we 
consider a bulk system with dephasing, turning to limited 
system size in the following section [6] 



5.1 Dephasing processes 

Perfect phase coherence is of course an idealization, and 
in real systems interference effects are killed by dephasing 
processes at some rate Jc/> = t7 . This means that the 
contribution of long diffusion loops to the interference- 
renormalized diffusion constant is suppressed. 

To incorporate dephasing processes in the self-consis- 
tent theory, we implement the simple prescription uj i— > 
uj + i/tip in (I39|) . This replacement can be shown to be ex- 
act within linear response theory for a rather large class 
of different systems with microscopic dephasing sources, 
such as electrons with spin-flip scattering by magnetic im- 
purites or photons being scattered by atoms with a degen- 
erate dipole transition '30]. Also mobile impurity atoms 
in a two-species experiment can be a source of dephasing 
randomness [31j . and in the limiting case of slowly mov- 
ing impurities an effective dephasing time is obtained that 
leads to the same qualitative consequences as the above 
substitution [32^. 

It is important to stress here that the intensity propa- 
gation kernel (PU)) keeps unchanged under dephasing, i.e., 
there u is not to be shifted by i/r^. This fact outlines an 
important difference between dephasing processes and ab- 
sorption. Dephasing kills interference, hereby affecting the 
self-consistent diffusion constant given by (l40|) . but con- 
serves particle number. Uniform absorption would cause 
particle losses at some constant rate 7^ = t~^ . It would 
not affect the relative phase between interfering ampli- 
tudes, but would damp their contribution as the loop size 
increases, thereby hindering the cumulative strength of 
interference effects. In the present paper, we consider con- 
servative dynamics and therefore do not need to discuss 
absorption, which could be conveniently incorporated in 
the theory by the substitution lu ^^ uj + iJTa both in (j20p 
and in ([55]) . 

In the present context of atoms in optical potentials, 
we consider that the main source of decoherence for the 
atom dynamics is the spontaneous emission of photons 
from the speckle field [8|9j . These spontaneous emission 
processes occur at a rate 



1 



vr 

m 



(49) 



where F is the natural linewidth of the excited atomic 
state and (5 is the detuning with respect to the atomic 
transition frequency. 

If the average time between successive spontaneous 
emission events is of the order of — or even shorter than — 
the Boltzmann transport time, t^ ;$ tb, then this intense 
dephasing will be shown below to result in purely classical 
transport already discussed in section [31 If on the contrary 
T,/, S> tb , then the regime of coherent diffusive transport 
is reached. Since t^/tb scales as l/|(5|/-r, this can be al- 
ways achieved at constant potential height provided the 
detuning |(5|/i^ is sufficiently large. 



5.2 Weak localization and residual diffusion 

After the replacement cj ^ cj + i/r^ in the self-consistent 
implicit equation (j42p . we can now discuss its solution 
along the lines of section 14.21 Let us consider the station- 
ary limit tj — > 0. The stationary diffusion constant solves 
now the self-consistent equation 



i?B(fc) 



1 



1 



TTfc^T 



■log 



1 



TBDBik) 



(50) 



In the case of intense dephasing r^ < tb, we get 
D -^ Db, a result easy to understand: coherent effects 
are scrambled at such a high rate that Boltzmann diffu- 
sion is never affected by interference. The dynamics thus 
stays entirely classical, and the results of section |3] apply. 
Of course, the random momentum recoil transferred to 
the atoms by rapid spontaneous emission of photons in- 
duces already an atomic diffusion by itself [33], but here 
we assume that this effect is negligible compared to the 
disorder-induced diffusion under study. 

More interesting is now the regime of phase-coherent 
transport r^ ^ tb- As in section 14.21 the short-time 
dynamics ujtb '^ 1 is again dominated by the classical 
Boltzmann diffusion. However, the limit tj — > does not 
lead anymore to strong localization in the strict sense. 
Instead, in the stationary limit, the residual dynamics 
stays diffusive. In the case of moderate dephasing, i.e. 
T^ ^ T-ioc = ^ioc/-^B, from ((50)) one finds the weak- 
localization correction to the diffusion constant in the form 



DwL(fc) = DBik) 



2TTm 



ln(l 



T0S 

tb' 



(51) 



In the last case of very weak dephasing, r^ ^ tioc , the 
residual diffusion constant is 



D{k) 



T0 



DB{k) 



Tloc(fc) 

T0 



(52) 



We can picture this regime by saying that coherent effects 
tend to localize the particle but phase-breaking mecha- 
nisms impose a residual slow diffusion with D <C Db'- 
the particle propagates coherently over a relatively long 
time and thus localizes in time rioc on a scale ^loc until a 
phase-breaking event occurs (typically every r^), resetting 
the diffusive Boltzmann dynamics followed by localization 
after another rioc, etc. 
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Fig. 4. Solid lines: diffusion constant D/Db given by eq. (|50|l 
as a function of k(^ for increasing disorder strength (from left 
to right) r] =0.05 (green), 0.1 (red), 0.2 (black). Decoherence is 
driven here by spontaneous emission induced by the light-shift 
speckle field itself, cf. 8 9 . The speckle detuning is chosen 
at 5 = 3.29 X lO^-T, where F is the width of the atomic ex- 
cited state, realizing highly phase-coherent propagation with 
t^/tb ~ 10^. The crossover from weak localization at high k 
to residual diffusion at small k is rather sharp. Dashed curves: 
corresponding plots of the ratio Dwh/Ds as given by (|5ip . The 
localization onset at fee where Dwl ~ gives a reasonably good 
approximation to the crossover. 



5.3 Localization onset in presence of dephasing 

Let us now examine more closely the crossover from weak 
to strong localization as a function of relevant parame- 
ters. Figure S] shows a plot of D/D-q solving ([50)1 (to- 
gether with the weak localization prediction ([5T|) ) as a 
function of kQ for different values of disorder strength rj. 
The laser detuning from the atomic resonance is held fixed 
at 5 = 3.29 X IQ'^F. For k(^ < 1, the low-energy limit ([28a|) 
then permits to estimate the effective dephasing rate as 
t,p/tb = 2Trri\d\/r, which varies slightly with rj from curve 
to curve but is exceedingly large (of the order of 10^) in 
all cases. 

Clearly, the curves show a sharp crossover between the 
weak localization regime at large momentum and resid- 
ual diffusion at small momentum. The crossover between 
these two regimes can be conveniently defined through the 
condition t^ = rioc. Recasting the weak-localization cor- 
rections (ISTI) in the form 
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WL 



1 



Db nkls 
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TB + Tloc 
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(53) 



one immediately sees that the crossover condition r^ = 
Tioc is reached at the point where the weakly localized 
diffusion constant would vanish: Dwl = 0. 

Solving for this condition defines an effective mobility 
edge fcc in momentum space as a function of the disorder 
parameters such as its strength r/ = VlEc^. Figure [5] shows 
a numerical plot of k^ as a function of t] as obtained for 
the 2D speckle correlator ([3]) in the case of spontaneous 
emission at detuning (5 = 3.29 x lO^-T. At low energy, one 
can use (I28al) to obtain: 
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Fig. 5. Critical wavenumber fccC (defining the effective onset 
of strong localization) as a function of the disorder strength 77 
in the optical speckle field with correlation ((3} and detuning 
& = 3.29 x VS' r . The dashed line is the low-energy behaviour 



This low-energy asymptotics is also shown in figure[5] Note 
that the weak-disorder condition 77 ^ fee C is well satisfied 
for r; <C 1. 

From a practical point of view, one can thus say that an 
effective strong localization regime is reached when only 
residual diffusion with D ^ D-q takes place. Consequently, 
one can earmark all momentum states with fc < fcc as 
"localized" and all states with fc > fee as "extended" . The 
localized fraction of atoms in this sense is 



Ploc 



fc<fc, 



dk 



^7ro(k)=/3^(2-/3^) 



(55) 



fccC = 2r; v/ln(l + 2^77|^|/r), [kC, « 1). 



(54) 



where /? = min(l, /cc/Zcmax)- Not surprisingly, localizing a 
large fraction of atoms requires that almost all momenta 
present in the initial distribution arc below the effective 
mobility edge. 



6 Finite-size systems 

The impact of long diffusion loops is also limited by spa- 
tial boundaries in finite-size systems. As justified in [34], a 
correct way to deal with finite-size systems is to introduce 
a space- dependent renormalized diffusion constant. A de- 
tailed calculation also requires to know which boundary 
conditions are imposed at the edges of the system. Obvi- 
ously, periodic or fixed boundary conditions (Dirichlet or 
Neumann) impose a discrete energy spectrum while open 
boundary conditions give a continuous spectrum and thus 
the impossibility of true localization (in the strict mathe- 
matical sense) at infinite time, because of leakage at the 
edges. 

A much simpler recipe is the replacement g^ — j g^ + 
1/L^ in the integrand of the return probability ((40l) (an 
equivalent alternative is to use 1/L as an infrared cut- 
off for the g-integral), thereby suppressing the contribu- 
tion of long loops. Without a precise prescription of what 
happens at the boundaries (absorption, reflection with a 
phase shift, periodic boundary conditions, etc.), this is 
the simplest approach. As we will see in the following, it 
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leads nevertheless to physically reasonable results in the 
localized regime. Of course, the system size must not be 
too small for diffusion to take place, and one should have 

In order to implement this simple approach, one may 
still work with the complex-frequency method used in the 
previous section and use the formal prescription to — *■ 
oj + i/tl with Ti = L^ /D, a form which is sugges- 
tive of the average time it would take an atom with ef- 
fective diffusion constant D[ijj,k) to reach the bound- 
ary. This cutoff time is obtained from the Thouless time 
TTh = L^/Db if one replaces the bare diffusion constant 
Db by the self-consistently renormalized D. Note that al- 
though tl = L^ / D formally acquires an imaginary part if 
D{lo, k) does, the replacement in ((40|) still provides a real 



l/L^ because D{u},k) 



momentum replacement q^ —> q'^ 
is divided out. 

In principle, one should also use the (non-self- 
consistcnt) replacement to —t lo + j/rxh in the intensity 
propagator pO|l. Like uniform absorption, this leakage 
would lead to a global exponential damping of the prob- 
ability density, also known as the fundamental Holstein 
mode of the system [35]. In the localized regime where 
the system size is larger than the bulk localization length, 
L 3> Cioc(fc), this damping (whose precise effect would 
again depend on the boundary conditions) is expected to 
be small and can be neglected. 



6.1 Onset of localization in a finite-size system 
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Fig. 6. Reduced diffusion constant D{k)/DB{k) as a func- 
tion of reduced atomic wave number fc^ for increasing disorder 
strength (from left to right) rj =0.05 (green), 0.1 (red), 0.2 
(black) . Dashed line: only the finite size L = 4 mm of the sys- 
tem is taken into account. The onset of localization with van- 
ishing D happens when the localization length of the infinite 
system equals the system size. Solid line: when dephasing ef- 
fects such as spontaneous emission (detuning 5 = 3.29 x 10^ F) 
are taken into account, the transition is smoothed. 



Note that this correction is algebraic in i/<^ioc while ab- 
sorption at the boundary is expected to be exponentially 
small. Les us repeat that this simple approach is valid 
only if the system size is larger than the bulk localization 
length. When this condition is violated, a space-dependent 
diffusion constant and proper boundary conditions must 
be considered. 



We now briefly discuss the case of a finite-size system 
in the absence of dephasing mechanisms (r^ -^ oo). Us- 
ing the finite-size replacement cu —^ to + i/TL in the self- 
consistent equation (|42p one determines the stationary dif- 
fusion constant in the limit uj ~f 0, which turns out to be 
a non-zero real quantity when ^ioc(fc) > L : 
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(56) 



As this is precisely the regime where boundary conditions 
must be carefully specified, this result must be taken with 
a grain of salt, and only gives an indication of how much 
the classical transport is slowed down by interference ef- 
fects. 

The situation is more interesting below the critical mo- 
mentum kc defined by the condition ^ioc(fcc) = L. There, 
the stationary diffusion constant vanishes, implying local- 
ization, and the asymptotic low-frequency behaviour is: 



D{u;,k)^-tuj^{l{k). 
with a modified localization length [9]: 

11 1 

~ l2- 



e(fc) ec(fc) 



(57) 



(58) 



This expression says that cutting off very long loops in 
a finite system slightly increases the localization length. 



6.2 Localization in a finite-size system with dephasing 



We account for the possibility of transport occuring in- 
side a finite domain while suffering from dephasing by 
summing the respective rates, Tq = t^ + r 7 , a pre- 
scription reminiscent of the Matthiessen rule [20] . If one 
prefers to reason in terms of length scales, one may de- 
fine the effective (complex) phase-coherence length ^^ = 
y^Dr^ and use L'^^ = [L~^ + ^7^]^^^ as momentum cut- 
off through q^ ^ q^ + L^^. Since the complex diffu- 
sion constant is divided out, note again that this pre- 
scription is just equivalent to use real cutoffs through 



Dq' 



r-'+D{q^ + L-^). 



If both finite size and dephasing are taken into account, 
there is no strict onset for localization, but, as in the sit- 
uation with pure dephasing, an effective threshold below 
which residual diffusion is tiny, D ^ Db ■ Figure [6] shows 
the ratio of the diffusion constant — obtained by solving 
the self-consistent equation ([1^ — to the Boltzmann diffu- 
sion constant, when the finite size of the system is taken 
into account, and when dephasing effects are included. 
Comparison with fig. 3] shows that the effect of the finite 
system size in this case is rather small. 
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Fig. 7. Solid curves: ratio of momentum-averaged diffusion 
constants V/Vb as a function of fcmaxC of sm expanding 
BEC with momentum distribution (I25p for increasing disor- 
der strength (from left to right) rj = 0.05 (green), 0.1 (red), 0.2 
(black). The speckle correlation function is given by Q, the 
detuning has been fixed to S = 3.29 x 10^ F. Dashed curves: 
INvl/Cb obtained from the weak- localization prescription (|51|l 
and (f5^ . 



7 Localization in experiments 

7.1 Effect of the initial momentum distribution 

Knowing the diffusion constant D{k) for each momentum 
k, the average diffusion constant (P2)) can be computed by 
integration over the Thomas-Fermi momentum distribu- 
tion ((25|) of an expanding BEC. One can use either the 
self-consistent diffusion constant D{k) or its approxima- 
tion in the weak localization regime D^i^(k). In the latter 
case, the integral must be cut below kc in order to avoid 
counting an unphysical negative D^f^i,: 
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fe>fcc 
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^7ro(k)i:)wL(fc)- 



(59) 



Fig. [7] shows the ratio V/V-q as a function of fc^axC foi' 
different values of 77 and the 2D speckle correlation ^ . A 
value of T>/T>B close to 1 means that quantum corrections 
to transport are difficult to observe. Conversely they be- 
come very strong when this value approaches zero. Atoms 
are then practically completely localized. The threshold in 
k to observe localization increases with r/, or equivalently, 
with kc- We find again the criterion fcmax ~ ^c to observe 
"strong localization" rather unambiguously. 

In [13], a Gaussian initial wavepacket is discussed in- 
stead of a Thomas-Fermi distribution, leading to a Gaus- 
sian momentum distribution centered at fc = with dis- 
persion Ak. Fig. [8] gives our results for the ratio V/Vb as 
a function of AkC, for different values of -q. Compared to 
Fig. [71 there are no qualitative changes and no big quan- 
titative changes. Hence the conclusions keep the same. 



7.2 Experimental realizability 

In order to estimate the orders of magnitude expected 
for a realistic 2D experiment, we will use the various 
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Fig. 8. Same as in figure [T] but for a Gaussian momentum 
distribution with dispersion Ak. 



experimental parameters of the quasi-lD Orsay experi- 
ment [5]. There, a BEC condensate of Rb*^ atoms {m k, 
1.4410~^^kg) is produced with a healing length ^ « 
0.40 /xm. The condensate evolves in an optical speckle with 
correlation length 1^ « 0.26 ^m such that fcmaxC = C/C ~ 
0.65. The optical speckle size L « 4 mm is far-detuned 
from the Rb transition frequency by 5/ F « 3.29 x 10^. 
The speckle fluctuation strength can be adjusted by vary- 
ing the laser power, up to 77 « 0.046. For reasons explained 
below, we performed calculations for r? = 0.05, 0.1 and 0.2. 
Since "q/kma^C is then at most 0.3 (and therefore suffi- 
ciently less than 1), we can use the weak-disorder results, 
keeping in mind that exact numerical values can differ 
from the ones computed in lowest-order perturbation the- 
ory. 

Compared to the ID situation, the main difficulty in 
2D is to reach a sufficiently short transport mean free path 
£b- Indeed, the exponential dependence of the localization 
length ^loc with Ib — see eq. ([44]) — implies that kis must 
not be much larger than unity, so that fc^ioc itself can be at 
most of the order of few hundreds or thousands, meaning 
^loc smaller or comparable to 1 mm. Then, the Boltzmann 
diffusion constant, eq. ([^5)l . will be of the order of a few 
h/m, that is a few 10~^m^s^^. As the lifetime of an ultra- 
cold atomic gas is typically limited to a few seconds, it 
is clear that the maximum diffusive expansion of the gas 
will be limited to a few hundred fiui. This in turn implies 
that localization by interference effects will be observed 
only if the localization length is decreased to a sufficiently 
small value, of the order of a fraction of a mm. This is 
very different from the ID situation, where localization 
can be observed even for a large transport mean free path 
(between 50 /im and 2 mm in [5]) and a large Boltzmann 
diffusion constant. A fortunate consequence is that the size 
of the speckle optical potential can be limited to a fraction 
of mm, making it possible to focus the laser power on a 
smaller area, such that values of 77 up to 0.2 should be 
reachable. 

With these parameter values, the effective mobility 
edge is predicted to be at fccC ~ 0.345, 0.571 and 0.818 
for 77 = 0.05, 0.1 and 0.2, respectively. We thus predict 
that about 48%, 95% and 100% of the atoms will be lo- 
calized, respectively. Evidence of the 2D localized regime 
should thus be easily observable. 
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With rj = 0.2, the maximum localization length (for 
atoms with k — fcmax) is about 48 /xm, and the localiza- 
tion time around 2.6 s. Loss of coherence by spontaneous 
emission shoud occur at a very small rate, of the order of 
t7^ « 6.7 X 10~^s~^, leading to a tiny residual diffusion 

constant, of the order of 10""' times the Boltzmann dif- 
fusion constant. These very favorable numbers are due to 
the choice of a far off-resonance laser beam to create the 
optical potential. A different technical choice for the laser 
beam could of course lead to a different situation. 

In conclusion, these estimates show that strong local- 
ization of matter waves in 2D speckle potentials is within 
experimental reach. 



8 Conclusion 



We have studied quantum diffusion of matter waves prop- 
agating in 2D optical speckle potentials within the frame- 
work of a self-consistent localization theory. We have con- 
sistently incorporated the effect of dephasing mechanisms 
and system boundaries. We have calculated the average 
probability density p(v,t) of matter waves propagating in 
a disordered potential for any initial Wigner phase-space 
distribution. This quantity allows to calculate the variance 
of the expanding cloud. We have shown that dephasing 
processes prevent true localization from taking place and 
induce a residual diffusion in the long-time limit. In certain 
regimes of parameters, the resulting stationary diffusion 
constant is so small that, for all practical purposes, one 
can define a localization onset by the momentum where 
the bulk weakly localized diffusion constant vanishes. This 
mobility edge then separates "localized" states from "ex- 
tended" ones, allowing for a practical definition of the lo- 
calized fraction of atoms. These quantities serve as exper- 
imental parameters to study weak and strong localization 
from a practical point of view. 

We have specialized our results to the case of an inter- 
action-driven BEC expanding in the presence of the disor- 
dered optical potential and to the case of a Gaussian mo- 
mentum distribution. In any case strong localization can 
only be achieved for cold enough atoms. Taking state-of- 
the-art figures we have shown that Anderson localization 
of matter waves in 2D optical speckle potentials should be 
observable in near-future experiments. 



The authors would like to thank Guillaume Labeyrie for raising 
this subject and for many helpful discussions about this work. 
This work was supported by the DA AD, the BFHZ-CCUFP 
and the Marie Curie program (contract number HPMT-2000- 
00102) at Laboratoire Kastler Brossel (Universite Pierre et 
Marie Curie / Ecole Normale Superieure, UMR 8552 du 
CNRS). RK would like to thank Professor B.-G. Englert and 
the Quantum Information Lab at the National University of 
Singapore for their kind hospitality (A*STAR Temasek Grant 
012-104-0040). 



References 

1. M. Lewenstein, A. Sanpera, V. Ahufinger, B. Damski, 
A. Sen, U. Sen, Adv. Phys. 56, 243 (2007) 

2. P.W. Anderson, Phys. Rev. Lett. 109, 1492 (1958) 

3. B.v. Tiggelen, in Diffuse Waves in Complex Media, edited 
by J.P. Fouque (Kluwer, 1999), pp. 1-60 

4. B. Kramer, A. MacKinnon, Rep. Progr. Phys. 56, 1469 
(1993) 

5. J. Billy, V. Josse, Z. Zuo, A. Bernard, B. Hambrecht, P. Lu- 
gan, D. Clement, L. Sanchez-Palencia, P. Bouyer, A. As- 
pect, Nature 453, 891 (2008) 

6. G. Roati, C. D'Errico, L. Fallani, M. Fattori, C. Fort, 
M. Zaccanti, G. Modugno, M. Modugno, M. Inguscio, Na- 
ture 453, 895 (2008) 

7. E. Abrahams, P.W. Anderson, D.C. LicciardcUo, T.V. Ra- 
makrishnan, Phys. Rev. Lett. 42 (1979) 

8. R.C. Kuhn, C. Miniatura, D. Delande, O. Sigwarth, C.A. 
Miiller, Phys. Rev. Lett. 95, 250403 (2005) 

9. R. Kuhn, O. Sigwarth, C. Miniatura, D. Delande, C.A. 
Miiller, New J. Phys. 9, 161 (2007) 

10. D. VoUhardt, P. Wolfle, Phys. Rev. B 22, 4666 (1980) 

11. B. Shapiro, Phys. Rev. B 25, 4266 (1982) 

12. L. Sanchez-Palencia, D. Clement, P. Lugan, P. Bouyer, 
G.V. Shlyapnikov, A. Aspect, Phys. Rev. Lett. 98, 210401 
(2007) 

13. B. Shapiro, Phys. Rev. Lett. 99, 060602 (2007) 

14. D. Clement, A.F. Varon, J. A. Retter, L. Sanchez-Palencia, 
A. Aspect, P. Bouyer, New J. Phys. 8, 165 (2006) 

15. J.W. Goodman, in Laser speckle and related phenomena, 
edited by J.C. Dainty (Springer- Verlag, 1975) 

16. M. Hartung, T. Wellens, C.A. Miiller, K. Richter, 
P. Schlagheck, Phys. Rev. Lett. 101, 020603 (2008) 

17. M. HiUery, R.F. O'Connell, M.O. Scully, E.P. Wigner, 
Phys. Rep. 106, 121 (1984) 

18. J. Rammer, Quantum Transport Theory (Perseus Books, 
Reading, Mass., 1998) 

19. Y. Kagan, E.L. Surkov, G.V. Shlyapnikov, Phys. Rev. A 
54, R1753 (1996) 

20. N. Ashcroft, D. Mermin, Solid State Physics (Saunders 
College, Philadelphia, 1976) 

21. E. Akkermans, G. Montambaux, Mesoscopic physics of 
electrons and photons (Cambridge University Press, Cam- 
bridge, 2007) 

22. S. Chandrasekhar, Radiative Transfer (Dover Publica- 
tions, New York, 1960) 

23. V. M. Apalkov, M. E. Raikh and B. Shapiro, J. Opt. Soc. 
Am. B 21, 132 (2004) 

24. E. Akkermans, G. Montambaux, J.L. Pichard, J. Zinn- 
Justin, Mesoscopic Quantum Physics (Proceedings of the 
Les Houches Summer School Session LXl) (Elsevier Sci- 
ence, Amsterdam, 1995) 

25. O.I. Lobkis, R.L. Weaver, Phys. Rev. E 71, 011112 (2005) 

26. L. Gor'kov, A. Larkin, D. Khmel'mtskii, Pis'ma Zh. Eksp. 
Teor. Fiz. 30, 248 (1979) [JETP Lett. 30, 228 (1979)] 

27. G. Bergmann, Phys. Rep. 107, 1 (1984) 

28. M. Abramowitz, I. Stegun, Handbook of mathematical 
functions with formulas, graphs, and mathematical table 
(Dover Publications, New York, 1972) 

29. S. E. Skipetrov, A. Minguzzi, B. A. van Tiggelen, B. 
Shapiro, Phys. Rev. Lett. 100, 165301 (2008) 

30. C.A. Miiller, C. Miniatura, E. Akkermans, G. Montam- 
baux, J. Phys. A: Math. Gen. 38, 7807 (2005) 



C. Miniatura et al.: Quantum Diffusion of Matter Waves in 2D Speckle Potentials 13 



31. S. Ospelkaus, C. Ospelkaus, O. Wille, M. Succo, P. Ernst, 
K. Sengstock, K. Bongs, Phys. Rev. Lett. 96, 180403 
(2006) 

32. E. Kogan, Eur. Phys. J. B 61, 181 (2008) 

33. G. Grynberg, C. Robilliard, Phys. Rep. 355, 335 (2001) 

34. N. Cherroret, S.E. Skipctrov, Phys. Rev. E 77, 046608 
(2008) 

35. G. Labeyrie, E. Vaujour, C.A. Miiller, D. Delande, 
C. Miniatura, D. Wilkowski, R. Kaiser Phys. Rev. Lett. 
91(22), 223904 (2003) 



